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ABSTRACT 

We study gravitational lensing by the family of scale-free galaxies with flat rotation 
curves. The models are defined by a shape function, which prescribes the radius of 
the isophote as a function of position angle from the major axis. The critical curves 
are analytic, while the caustic network is reducible to a simple quadrature. The cusps 
are always located at the turning points of the shape function. We show that the 
models with exactly elliptic isophotes never admit butterfly or swallowtail cusps and 
so have at most 4 (or 5) images. Higher order imaging is brought about by deviations 
of the isophotes from pure ellipses - such as pointedness caused by embedded disks or 
boxiness caused by recent merging. The criteria for the onset of sextuple and octuple 
imaging can be calculated analytically in terms of the ellipticity e and the fourth-order 
Fourier coefficients (04 and 64) used by observers to parametrise the isophote shapes. 
The 6 or 8 images are arranged roughly in a circle, which appears as an incomplete 
Einstein ring if inadequately resolved. Using data on the shapes of elliptical galaxies 
and merger remnants, we estimate that ~ 1% of all multiply imaged quasars may be 
sextuplet systems or higher. Forthcoming satellites like the Global Astrometric Inter- 
ferometer for Astrophysics (GAIA) will provide datasets of ^ 4000 multiply imaged 
systems, and so 40 will show sextuple imaging or higher. 
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1 INTRODUCTION 

Almost all the known sixty or so gravitational lens sys- 
tems are 2 or 4 image configurations (see Pospieszalka et 
al. 1999 for details of the gravitational lensing database 
which maintains a list of candidates) . There are two systems 
in which the presence of a weak, central image has been 
claimed, namely APM08279+5255 (Ibata et al. 1999) and 
MG1131-I-0456 (Chen & Hewitt 1993), but these detections 
are controversial. There are only two lens systems which are 
definitely established to have more than 4 images of a single 
source. The first is CL0024-I-1654 which comprises 8 images 
of a single blue galaxy lensed by a foreground cluster of 
galaxies (CoUey, Tyson & Turner 1996; Tyson, Kochanski & 
Deir Antonio 1998). The second is B1359-|~154, which con- 
sists of 6 images of a single radio source (Rusin et al 2001). 
Here, the lens is believed to be a compact group of galaxies. 

For non-singular lenses, it is well-known that the total 
number of images is odd and that the number of even parity 
images exceeds the number of odd parity images by one 
(e.g., Burke 1981; Schneider, Ehlers & Falco 1992, chap. 5). 
When the source lies within the central or tangential caustic. 



such lenses typically give rise to 5 images. Four of these lie 
roughly on an Einstein ring, while the fifth is highly de- 
magnified and located near the centre of the lensing galaxy. 
For cusped lenses, the convergence (or surface density in 
units of the critical density) k behaves Uke 



lim K{r) = 0(r 

r— >0 



(1) 



and so diverges at the galaxy centre. Evans & Wilkinson 
(1998) showed that there are 4 images if the central cusp is 
stronger than isothermal (7 > 1), and 5 images if the cusp is 
weaker than isothermal (7 < 1), again assuming the source 
lies within the central caustic. Hence, the observed abun- 
dance of 4 image systems may be caused either by strong 
cusps (in which case there are only 4 images) or by small 
core radii (in which case the fifth image is so de-magnified 
as to be invisible against the much brighter lensing galaxy) . 

The purpose of this paper is to ask whether sextuplet 
or octuplet images may exist? If so, what properties of the 
lensing galaxy cause such higher order imaging? Note that in 
this paper, by a slight abuse of terminology, we use sextuple 
imaging to refer to systems with 6 or 7 images; a central 
seventh image may exist if the galaxy is soft-centered but 
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will in general not be discernible. Similarly, octuple imaging 
may refer to 8 or 9 images. Systematic programs like the 
Cosmic Lens All-Sky Survey or CLASS (e.g., Myers et al 
1995, Myers 1999) are now providing a much larger dataset 
of multiply-imaged quasars than the earlier serendipitous 
discoveries. The future holds out even brighter prospects, 
as the Next Generation Space Telescope (NGST) and the 
Global Astrometric Interferometer for Astrophysics (GAIA) 
will discover many thousands of such systems (e.g., Barkana 
& Loeb 2000; de Boer et al 2000). Hence, there are golden 
opportunities for uncovering comparatively unusual image 
configurations. 

Section 2 provides an introduction to the lensing prop- 
erties of scale-free galaxies with flat rotation curves, n We 
show that the caustic network, together with the locations 
of the cusps, can be worked out exactly. Section 3 considers 
two simple and very familiar models, namely those with el- 
liptic isophotes and elliptic equipotentials (e.g., Kassiola & 
Kovner 1993). Such models, if acting as isolated lenses, are 
unable to produce higher order imaging. Section 4 identifles 
the physical causes of higher order imaging in terms of de- 
viations of the isophotal shapes of the lensing galaxy from 
pure ellipses. Simple models are provided for lenses with 
boxy, disky and squarish isophotes and the onsets of sextu- 
ple and octuple imaging are calculated analytically. Finally, 
Section 5 concludes with a discussion of the frequency with 
which such higher order image systems occur, together with 
an assessment of forthcoming observational opportunities. 



2 SCALE-FREE MODELS WITH FLAT 
ROTATION CURVES 

2.1 Mass and Potential 

Scale- free galaxy models with flat rotation curves are widely 
used in galactic astronomy and dynamics (e.g., Toomre 1982; 
Evans 1993; Evans, CaroUo & de Zeeuw 2000). The isophotes 
of a scale-free galaxy have the same shape at every radius. 
The isophotes are completely described by a shape function 
G{9), which depends only on the position angle 9 with re- 
spect to the major axis. In such galaxies, the convergence k 
and deflection potential <j!> satisfy 



2r ' 



= rF{e), 



(2) 



where (r, 0) are familiar polar coordinates in the lens plane. 
The surface density falls off along any ray like the reciprocal 
of distance and so the three-dimensional density falls off like 
the inverse square of distance. The models are projections 
of axisymmetric and triaxial generalisations of the familiar 

* Such models are often called "isothermal" in the lensing litera- 
ture. This is an unfortunate nomenclature, as the velocity distri- 
butions are not in general isothermal. We prefer to use the more 
cumbersome term "scale free galaxies with flat rotation curves" 
to describe this family of models. The circular velocity curve of 
the spherical and axisymmetric models is everywhere flat. For the 
triaxial models, the circular velocity curve is not strictly speaking 
defined. Nonetheless, the structure of velocity space (and hence 
the elliptic closed orbit in the principal plane on which any cold 
gas lies) is the same independent of position. Hence, these models 
too have everywhere flat rotation curves. 



isothermal sphere. As \7'^(j> = 2k (see e.g., Schneider et al. 
1992), it is straightforward to establish that 



G{e) = F{e)+F"{e). 



(3) 



Given an arbitrary shape function G{9), it is possible to 
solve (^ by the method of variation of the parameters (e.g., 
Bronshtein & Semendyayev 1998, section 3.3.1.3.4) to obtain 

F{e)=sme G(i9)cos-!9di?-cose / G(i9) sin i9 di9. (4) 

The deflection angle has components 

G{-d)sm-dd^, (j}y= I G(i?)cosi?di9. (5) 
Jo 

This equation relates the lensing properties of the galaxy 
directly to the shape of the isophotes G(6). 



2.2 Critical Curves and Caustics 

For scale-free models with flat rotation curves, the magnifi- 
cation is simply related to the convergence via (e.g., Witt, 
Mao & Keeton 2000) 



det J = 1 - 2k. 



(6) 



The critical curves are given by the vanishing of the Jacobian 
and so they have the polar equation 



r = G{e) = F{e) + F"{e). 



(7) 



Let the Cartesian coordinates of the caustic be (C,??)- Then, 
the caustic is given through the lens equation as 

^ = F" (9) cos 9 + F' (9) sin 9, 

ri = F" (9) sin 9 - F' (9) cos 9. (8) 

If we introduce polar coordinates (s, a) in the source plane, 
then the equations (^) become 

scos{9 ~a) = F"{9), s sin{9 ~ a) = F' (9) . (9) 
By substituting (^) into we obtain the result 

C = G{9) cos 9+ [ G(i?)sini9di9, 



G(6l)sin6i 



G(i9)cosi9dj9. 



(10) 



Although this equation is simple to derive, it is rather re- 
markable. For scale-free galaxies with flat rotation curves, 
the caustics can be computed by a simple quadrature for any 
isophotal shape. Generally speaking, in gravitational lensing 
it is very hard to establish the maximum number of im- 
ages that can be generated by a given lens. Often, numeri- 
cal searches are the only way of tackling such questions and 
they risk missing small but important domains of parameter 
space. By contrast, here the caustic network is directly avail- 
able to us and the onset of sextuple and octuple imaging can 
be found exactly. 



2.3 The Role of Cusps 

The appearance of swallowtails or butterflies usually leads 
to areas in the source plane where 6 or more images are 
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possible. Knowledge about the position of cusps is therefore 
very important. Since the caustic is given by the parametric 
representation (^), it is straightforward to derive the cusp 
locations. At a cusp, the caustic stops and reverses direction, 
so that the derivative of the caustic coordinates with respect 
to the parameter must vanish: 



de 



G'{e)cose = 0, 



G'{e) sine = 0. 



(11) 



Both equation can only vanish simultaneously if the common 
factor G'{9) vanishes. Therefore, the location of the cusps 
^cusp is simply given by the equation G'(Scusp) ~ 0. This fact 
is remarkable since usually the equation for the cusps is very 
complicated. The cusps can occur only at the turning points 
of the shape function, so they are related to the properties 
of the isophotes in a straightforward way for these models. 

In models without higher order cusps - such as swal- 
lowtails and butterflies - the caustic only has 4 simple cusps 
or folds. If a galaxy model permits the occurrence of swal- 
lowtails or butterflies, then it can show sextuple or octu- 
ple imaging. The onset of swallowtails happens when the 
second derivative of the shape function vanishes as well, 
G"(6lcusp) = G"'(6lcusp) = (see Schneider et al. 1992; Kee- 
ton, Mao & Witt 2000) so the flrst non-zero derivative is 
the third. Similarly, the onset of butterflies happens when 
all of the first three derivatives of the shape function vanish, 
G"(6lcuBp) = G"(6icusp) = G"'(6icusp) = 0, so the first non-zero 
derivative is the fourth. 

The zeros for the cusps Scusp are related to the length 
of the caustic. For the length of the caustic, we may write 



VWWf + WWVde = / \G'{e)\de 



lea 



"cUBp/2 

= 2 ^ ^ \G(6cusp.2i) ^ G{9cusp,2i-l)\- (12) 
1=1 

Finally, the area under the caustic can also be given 

c2tv 

Acau.tic = / i{e)r^'(e)de 



de I di;G'{e)G'{ip) cose sin 4i. (13) 
Jo 

This holds good provided the caustic does not self-intersect. 



2.4 Area and Mass 

For the spherical case (G(0) — const), it is well known that 
the mass inside the Einstein ring is directly related to the 
Einstein radius by Me ~ r| (cf. Schneider et al. 1992). 
The area enclosed by the Einstein ring (i.e. the area inside 
the critical curve) is related to the mass by tt Merit. curve = 
^crit. curve. From au observational perspective, the relation 
becomes tt Merit. curve = Seritl^dArit. curve, whcrc the area is 
measured in arcsec^ and the critical surface mass density is 
given by 



(14) 



47rGi3dI>ds 

Here, D^,Db and Dds are the distance to the deflector, the 
distance to the source and the distance from the deflector 




Figure 1. Critical curve (solid line) and reciprocal critical curve 
[as defined in eq. jl^)] (dashed line) for an elliptical density dis- 
tribution [q = 0.8) is shown. In both cases the mass inside the 
curve is proportional to the area, i.e. A = ttM. The filled circles 
show the positions of the images when the source is located at 
the origin. 




Figure 2. The isophotes of the family of models with shape func- 
tion ( ^ ) with q = 0.75. Note that a = 0.5 corresponds to ellip- 
tical isophotes, while a = 1.5 corresponds to elliptical equipoten- 
tials. All these lenses are unable to show higher order imaging. 



to the source respectively. Given the redshift of the lens and 
the source, this (|l^ ) enables the projected mass within the 
Einstein ring to be calculated for a spherical lens. In general, 
these useful relations do not remain valid for arbitrary, non- 
spherical lenses. 

However, we now show that they do remain true for 
any scale-free potential obeying eq. (□). The mass within 
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Figure 3. The onset of sextuple imaging as a function of the ellipticity of the isophotes e and the parameter a4 used to measure boxy 
and disky deviations. Upper panel additionally shows data from Saglia et al. (1993); triangles are elliptical galaxies, crosses are SOs and 
squares are SBOs in the Coma cluster. Lower panel additionally shows data from simulations of merger remnants by Heyl et al. (1994). 



the critical curve is 



Merit, cur' 



1 

1 

2^ 



27r l-RiO) 



K.{r, 9)rdrd9 



-'0 

2w 



G\e)de. 



The area inside the critical curve is 



A 



crit. curve — 



2-K 



R(0) 



rdrdO = 



G^{e)de, 



(15) 



(16) 



where we have used eq. (^ for the dependency of the critical 
curve on the angle. So, we obtain again vr Merit ~ ^crit- 

Many authors have suggested that the mass inside the 
ellipsoidal ring formed by the images of a quadruple lens 
is approximately constant, almost independent of the lens 
model; for example, the models of Kochanek (1991), Rix, 
Schneider & Bahcall (1994), Witt, Mao & Schechter (1995), 
Chae, Turnshek, & Khersonsky (1998) and Hunter & Evans 
(2001) for the Einstein Cross (Q2237+0305) are very differ- 
ent and predict different flux ratios for the 4 images, but 
they all agree on the projected mass within the central 0.9" 
enclosed by the images. We want to shed some light on this 
issue! Let us define the reciprocal critical curve as 



R{e) = c/rcrit.cur 



c/G{e). 



(17) 



The mass inside the reciprocal critical curve is 



M 



rccip. curve 



R(0) 



n{r, 9)rdrd6 — c. 



(18) 



and it too must remain constant. Let us choose the constant 
c in such a way that the area inside the reciprocal critical 



curve A 
ing 



rccip. curve 



2tt 



G~^{e)de 



/tt — CIS proportional to the mass, obtain- 



(19) 



As shown in Figure ^ we can imagine that the images are 
confined by the critical curve and the reciprocal critical 
curve. If the source position is closer to the caustic, then 
the images are located closer to the critical curve. In other 
cases, the images may be closer to the reciprocal critical 
curve. The image positions follow (approximately) the criti- 
cal curve or the reciprocal critical curve. The area within the 
critical and inverse critical curves is directly proportional to 
enclosed mass. Hence, to a good approximation, the area en- 
closed by the four images is as well. (Small deviations might 
arise due to the interpolation curve R(0) chosen to connect 
the images of the quadruplet). 



© 2001 RAS, MNRAS 000, 000-000 



Are there sextuplet and octuplet image systems? 5 



0.05 







-0.05 - 




(a)- 



-0.05 



0.1 
0.05 



-0.05 
-0.1 



0.05 



(c)- 






-0.1 -0.05 0.05 0.1 



Figure 4. The caustic network and isophote shape for elliptical 
galaxy models with (a,b) boxy and {c,d) disky isophotes. Butter- 
fly caustics appear on respectively the major and minor axes pro- 
jected onto the source plane. A pure elliptical isophote is shown 

for comparison in dotted lines, so that the boxy and disky isopho- 
tal deviations can be discerned. [Panels (a,b) are for a model 
with e = 0.05 and 04 = —0.02, panels (c,d) are for e = 0.1 and 
04 = 0.03]. 
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Figure 5. As the distortion increases, the butterfly cusps merge. 
If the source lies in the region where the two caustics overlap, 
then eightfold imaging occurs. [This shows a model with e = 0.1 
and 04 = -0.075]. 



3 SIMPLE MODELS 

3.1 Elliptic Isophotes and Equipotentials 

There are two very familiar models that fall into the class 
of scale-free galaxies with flat rotation curves. The first has 



elliptic isophotes, so that the shape function is 

G{9) oc (i(l + g-') + i(l-g-')cos26l)"' 

oc (cos^ 6 + q^^ sin^ 9^ ' . 
It is easy to find the deflection potential as 
(/> = a;0x + 2/0!/, 
with 



(20) 



(21) 



G(i?) sin-dM 



tan 



l-g2 



g2 cos2 e+sin^ 6 



cos 6 



G(t?)cost?dt9 



tanh 



l-g2 



q-^ COS'' ^ + sin 



(22) 



This result has been given before (see e.g., Kassiola & 
Kovner 1993; Kormann, Schneider, Bartelmann 1994). By 
explicit construction, Keeton et al. (2000) showed that an 
elliptical density distribution plus external shear can pro- 
vide 6 or 8 images configurations. This happens when the 
external shear is strong and aligned along the minor axis of 
the lensing galaxy. 

The second has elliptic oquijjotoutials. so that 



m « ^i(l + g-2) + i(l-9-2) 



cos 261 



' + sin 



(23) 



This model was introduced into lensing by Blandford & 
Kochanek (1987) and subsequently studied by others (e.g., 
Witt 1996; Witt & Mao 1997, 2000; Hunter & Evans 2001). 
The three dimensional analogues of these models arc well- 
known in galactic dynamics as power-law galaxies (Evans 
1993, 1994). It is straightforward to find the shape function 
as 



G{e)^q-^{\{l + q- 



^)cos 26)' 



(24) 



On geometric grounds, it seems likely that an elliptically 
stratified potential plus external shear cannot have more 
than 4 (or 5) images. The reason for this is the confinement 
of the image positions. As found by Witt (1996) and Witt 
& Mao (1997), the 4 (or 5) images of an elliptical potential 
are constrained to lie on a hyperbola. If any new image pairs 
appear, they must also lie on this hyperbola. Physically, we 
expect that each new image pair is located close to the Ein- 
stein ring. Since a hyperbola can have at most only 4 inter- 
section points with a ring-like curve, it is intuitively clear 
that an elliptical potential plus shear may not have more 
than 4 images (or 5 images if there is additionally one close 
to the centre of the lensing galaxy). An analytical proof of 
the maximum number of images for power-law galaxies with 
elliptic equipotentials is provided in Appendix A of Hunter 
& Evans (2001). 
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Figure 6. The onset of octuple imaging as a function of the el- 
lipticity of the isophotes e and the parameter 04 used to measure 
boxy and disky deviations. Data from simulations of merger rem- 
nants by Heyl et al. (1994) are overplotted. 



3.2 Number of Images 

Both these simple models are members of a more general 
class that have a shape function G(Q) of the form: 

G(6l)cx (|(l-f q-^)-f i(l-q-^)cos26l)"", (25) 

where a is constant. Figure |^ shows isophotes for three mem- 
bers of the family, namely a = 1/2 (elliptic isophotes), a = \ 
and a = 3/2 (elliptic equipotentials) . For the same g, the 
models become flatter with increasing a. The isophotes can 
become "peanut-shaped" ; this is valuable for modelling the 
peanut-shaped bulges of spiral galaxies (see e.g., Shaw 1987, 
Bureau & Freeman 1999). 

The positions of the cusps are found by solving G'{9) = 
and so 



mv 



n = 0, 1,2,3. 
It is straightforward to establish 



G"ie cusp 



)CX(1-Q-^) [l+g-2 + (_l)"(l. 



(26) 



(27) 



For flattened models {q 7^ 1), this never vanishes and so 
higher order cusps like swallowtails never occur. This pro- 
vides a simple proof that scale-free galaxies with flat rotation 
curves and elliptic isophotes or equipotentials have at most 
4 (or 5) images depending on the density singularity at the 
centre. Of course, the proof is slightly more general, holding 
good for all models with shape functions of form eq (Pq). 



4 ISOPHOTAL DEVIATIONS 

The isophotes of elliptical galaxies are well approximated by 
pure ellipses. Observers typically measure deviations from 
ellipses of the order of a few percent. This might betray 
the presence of an embedded disk or a central bar. Boxy or 
irregular isophotes may be caused by recent merging and 
interaction. Typically the difference between the isophote 
ri and a reference ellipse is expanded as a Fourier series 
(e.g., Binney & Merrifield 1998) 



y ^ a„ cos nd + bn sin n9, 



(28) 



n=l 



Provided the reference ellipse has been well-chosen, the 
Fourier coefficients are all small. The largest non-vanishing 
components are usually 04 (which is a measure of boxiness 
and diskiness) and 64 (which is a measure of skewness and 
squareness) . The 03 and 63 coefficients are significant only in 
elliptical galaxies with dust lanes. Observers normaUse the 
Fourier coefficients to the semimajor axis of the reference 
ellipse. Henceforth, we always use normalised Fourier coef- 
ficients, so that our models are directly comparable to the 
data. 

4.1 Boxiness and Diskiness 

Boxiness and diskiness has been of interest ever since Ben- 
der et al. (1989) pointed out possible correlations of the 
04 parameter with X-ray, radio and kinematic properties. 
This correlation remains poorly understood, especially as 
the boxiness and diskiness is partly a function of the incli- 
nation angle of the galaxy (Franx 1988). The fiattest galaxies 
have the largest values for |a4|. In fact, galaxies with an el- 
lipticity greater than 0.35 almost always have |a4| > 0.005 
(Bender et al. 1989). Approximately one third of elliptical 
galaxies show boxiness at the level a4 < —0.005, one third 
show diskiness at the level 04 > 0.005 and one third show 
deviations of \a4\ < 0.005. Values of |a4| much in excess of 
0.05 are uncommon for normal ellipticals. 
We study the shape function is 

G{e) oc [i(l + q"^) + i(l -q"^)cos26» -ylcos46i)]"^ .(29) 

When A vanishes, the model reduces to (|2^). Otherwise, 
if j4 > 0, then the model has disky isophotes, while if 
A < then the model has boxy isophotes. Observers typ- 
ically record the isophote shape in terms of the ellipticity 
e = 1 — q and the Fourier coefficient 04. In the limit of small 
deviations, A is related to the observables through the equa- 
tion 

A 



04 



32 



(4-3e)" + 0(e^A"). 



(30) 



These models can exhibit butterfly cusps at surprisingly 
small deviations. The condition for butterflies is G'{6^usp) = 



G"(0c 



G"'(5cusp) = 0, so that the first non-vanishing 



derivative is the fourth. This is satisfied at or beyond the 
critical value 



|j4crit,( 



(31) 



This condition can be re-written in terms of the observa- 
tional parameters 04 and e as 

1 e(4-3e)^(2-e) 



£14 



> 



256 



(1 



(32) 



This marks the onset of sextuple imaging. 

In Figure ^, the boundary between quadruple and sex- 
tuple imaging is shown in solid lines. In the upper panel, we 
also plot the data of Saglia, Bender & Dressier (1993) on a 
total of 58 elliptical galaxies, SOs and SBOs in the Coma clus- 
ter. In the lower panel, we plot the data on isophotal shapes 
of 424 merger remnants from simulations by Heyl, Hernquist 
& Spergel (1994). There are a number of things to notice. 
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Figure 7. The topology of the Fermat surfaces for instances of (left panel) sextuple and (right panel) octuple imaging. The even parity 
images are located at the minima of the Fermat surface, the odd parity images at the saddle points. Image locations are marked with 
filled circles. [The left panel corresponds to the model shown in Figure ^a) with source position ^ = 0.02, ?y = 0.0. The right panel 
corresponds to the model shown in Figure ^(a) with source position ^ = 0.005, »7 = 0.0]. 



First, very small values of |a4| suffice to produce sextuple 
imaging, if the galaxy is quite round. However, nearly round 
galaxies tend to have very small deviations from the pure el- 
liptical form. Second, all the elliptical galaxies in Saglia et 
al.'s (1993) sample lie in the quadruple imaging domain. 
The possibilities for sextuple imaging improve for SO and 
especially SBO galaxies. Saglia et al's sample contains only 
galaxies for which a stable value of 04 can be derived from 
the photometry and does not include galaxies which exhibit 
isophote twisting or those which show boxiness and diski- 
ness in the same isophotal contours. So, it probably under- 
estimates the frequency with which substantial deviations 
from pure ellipses occur. Nonetheless, out of a total of 58 
galaxies, 3 have the potential for sextuple imaging. Third, 
the merger remnants of Heyl et al. (1994) are typically boxy. 
This though is a rough rule-of-thumb, as some can appear 
disky. We see that 10 out of the total of 424 simulated merger 
remnants are irregular enough to permit sextuple imaging 
or higher. Bearing in mind that most of the lensing optical 
depth is in early-type galaxies, we reckon that the fraction 
of lenses that have the ability to cause sextuple imaging or 
higher is ~ 3%. 

Figure ^shows the caustic networks for two models. The 
upper panels refer to a model which has ellipticity e — 0.05 
and boxy isophotes with a4 — —0.02. The lower panels refer 
to a model which has e = 0.1 and disky isophotes with a4 = 
0.03. The locations of the cusps are given by the equations 



sin 2ecusp = and cos 20cusp = 



(33) 



This means that the 4 cusps remain on the major and minor 
axes with increasing distortion A, but 4 additional cusps 
appear as soon as the critical value Acrit.e is exceeded. Note 
the butterfly cusps develop on the major axis in the boxy 
case, and on the minor axis in the disky case. If the source 
lies within the butterfly caustic, then it is lensed into six 
images. 

The butterfly cusps touch at the onset of octuple imag- 
ing. Figure |5| shows an example when the butterfly cusps 
have merged. A source in the overlapping region will be 
lensed into eight images. The isophote shape is now very 
rectangular, but examples of such galaxies are known (see, 
for example, NGC 128 as pictured in Burbidge & Burbidge 
1959). The critical value of the distortion Acrit.g required for 
eightfold imaging is given by solution of the implicit equa- 
tion 



di9 cos i9 



[l(l+g-2) + i(l-g-2) cos2T?-lAcrit,8| cos4i?]i 



(34) 



As the distortion is always small, we can linearise this equa- 
tion and solve it to obtain 



2(/i - q) 



(35) 



This marks the onset of octuple imaging. The integrals I\ 



© 2001 RAS, MNRAS 000, 000-000 



8 N. Wyn Evans and Hans J. Witt 



and I2 are 



cos 



[|(l + g-2) + i(l-g-2)cos2,9]2 



itanh ^ -^/l — g2, 



0.05 







cos 1? COS 4i?di9 



/o [|(l + g-2) + i(l-g-2)cos2^]S 
= g[g" + 10gVl-4g(2 + g')Ji]/(l-g')^ (36) 

They can be evaluated numerically easily enough. Figure ^ 
shows the distortion required for octuple imaging in the 
plane of the observables. Overplotted is the simulation data 
of Heyl et al. (1994), from which it is clear that 5 of the 
merger remnants are sufficiently distorted to permit even 
eightfold imaging. 

Figure ^ shows contours of the Fermat potential 



(37) 



(38) 



G(e)« [i(l+g- 



(1 



) cos 261 - B sin 461)] ^ , (39) 
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for typical 6 image and 8 image configurations. Let us recall 
that the Fermat potential is proportional to the time delay 
(e.g., Schneider et al. 1992). The left panel shows the model 
of Figure ^a) with the source placed within the butterfly 
cusp so that it has 6 images (3 of even parity corresponding 
to the minima of the Fermat surface, 3 of odd parity corre- 
sponding to the saddle points). The right panel shows the 
model of Figure |^(a) with the source placed in the region 
of overlap of the butterfly cusps, so that there are now 8 
images, 4 of even parity and 4 of odd parity. Notice that all 
the images are at roughly the same distance from the centre 
of the lensing galaxy. In fact, the arrangement has obvious 
similarities with the deep Hubble Space Telescope images of 
CL0024-I-1654, which comprises 8 images of a single blue 
galaxy lensed by a foreground cluster (Tyson, Kochanski & 
Deir Antonio 1998), as well as the incomplete optical Ein- 
stein ring 0047-2808, which is the image of a high redshift 
early-type star-forming galaxy (Warren et al. 1999). 

Scale-free models with flat rotations have another ele- 
gant and unusual property. The relative time delay between 
the ith and jth image is simply related to the radial posi- 
tions ri and Vj via (Witt et al. 2000; see also Koopmans, de 
Bruyn & Jackson 1998; Zhao & Pronk 2001) 
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where Zd is the redshift of the deflector. Sextuple and octuple 
image systems typically have all the images in a ring-like 
configuration and so the relative time delays between any of 
the pairs is never very large. 

4.2 Skewed Squarishness 

Photometric studies of elliptical galaxies find that the other 
Fourier coefficient that is typically non-zero is hi (e.g., 
Franx, lUingworth & Heckman 1989). To study its effect, 
we examine the shape function 



where without loss of generality _B > 0. If the distortion 
B is non-vanishing, then the isophotes become skewed and 
squarish. This may be useful for modelling face-on bars in 
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Figure 8. The caustic network and isophote shape for galaxies 
with the shape function (^). Again, pure elliptical isophotes are 
shown in dotted lines for comparison. Both models have swal- 
lowtail cusps. [Panels (a,b) are for a model with e = 0.05 and 
64 = 0.02, panels (c,d) are for e = 0.05 and 64 = 0.06]. 
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Figure 9. The onset of sextuple imaging as a function of the 
ellipticity of the isophotes e and the parameter 64 used to measure 
skewed squarishness. 



spiral disks, rather like the lensing galaxy for Q2237+0305 
(the Einstein Cross). The observables are the ellipticity e = 
\ — q and the Fourier coefficient 64, which is 



64 = -(4- 
32^ 



3e 



''+o(e^ 



(40) 



Unlike 04, the sign of 64 carries no physical significance 
and can be inverted by a simple reflection. Observationally 
speaking, the 64 coefficient is usually, but not always, smaller 



© 2001 RAS, MNRAS 000, 000-000 



Are there sextuplet and octuplet image systems? 9 



than 04 (see e.g., Franx et al. 1989). The 64 coefficient has 
received less attention than 04, as it does not appear to 
be correlated with radio or X-ray properties. There are no 
extensive tables of data on 64 available in the literature. 
Nonetheless, individual galaxies with 64 as large as 0.05 are 
known, such as NGC 1700 (Franx et al. 1989). 

The condition for swallowtails is G'(^cusp) = 
G"(Scusp) = 0, so that the first non-vanishing derivative is 
the third. This is satisfied whenever B exceeds the critical 
value given 



-Bcrit.e = =1=^(1 - q 



(41) 



Again, it is useful to rewrite this in terms of the observables 
as 

1 e(4-3e)^(2-e) 



&4 



128 (1 - e)2 



(42) 



This marks the onset of sextuple imaging. 

Figure ^ shows the caustic network for a case when the 
swallowtails have just developed (e = 0.05 and 64 = 0.02), 
together with a case in which the swallowtails have over- 
lapped to create a domain in the source plane where eight- 
fold imaging occurs (e — 0.05 and 64 = 0.06). The cusps 
occur at the locations 



sin 261 



1 



(1 



256B 



2)2 I 
+ -. 



(43) 



The cusps move further away from the major and minor 
axes for increasing distortion B. We obtain 4 additional 
cusps when the swallowtails form as B reaches the crit- 
ical value -Bcrit,6. This happens at the angular locations 
Scusp = f (2n -I- 1), where n is an integer. If the distortion 
B increases still further, then the swallowtails can overlap. 
This is when eightfold imaging becomes possible. The crit- 
ical value i3crit,8 can be determined by requiring that the 
two inner pair of cusps just touch each other. Let us la- 
bel the cups in Figure M in numerical order counterclock- 
wise starting with 1 on the extreme right hand side. Then, 
the condition can be expressed as 5(&cusp,2) ~ C(^cusp,7) 
and r;(^cusp,2) = 7)(6'cusp,7). Simultaneously, the 6th and 8th 
cusps touch each other as well. The onset of octuple imaging 
can be calculated analytically, but it is rather cumbersome 
and so we resort to the computer. 

Figure |^ shows the domains in which sextuple and oc- 
tuple imaging is possible as a function of the ellipticity e and 
the Fourier coefficient 64. Unfortunately, tables of data of 64 
are not available, although 64 is typically smaller than 04 
and rarely exceeds 0.05. On comparing Figure ^ to the ear- 
lier Figures ^ and ^ we see that the skewed squarishness is 
less efficient at causing higher order imaging than diskiness 
or boxiness. 



5 CONCLUSIONS 

Many sextuple or octuple images of multiply lensed quasars 
are likely to exist. If a galaxy with a flat rotation curve has 
exactly elliptic isophotes, then it yields at most 4 or 5 images 
(depending on the behaviour of the convergence or surface 
density at the centre). Higher order imaging is caused by 
deviations of the isophotes from pure elliptic form, such as 



boxiness or diskiness. Even isolated boxy or disky ellipti- 
cal galaxies provide caustic networks with swallowtail and 
butterfly cusps and so permit sixfold or eightfold imaging 
if the source is located propitiously. The rounder the lens- 
ing galaxy, then the smaller the distortion of the isophotes 
required. The criteria for the onset of sextuple and octuple 
imaging can be given exactly in terms of the ellipticity e 
and the Fourier coefficients 04 and 64 used by observers to 
measure deviations such as boxiness, diskiness and skewed 
squarishness (see eqs. ^ |35[ ^). 

The first examples of galaxy models that can produce 
sextuple or octuple images were found by Witt & Mao (2000) 
and Keeton et al. (2000). These models have elliptical den- 
sity contours with the external shear perpendicular to the 
flattening. These are rare cases, as they require large values 
of the external shear for the cross-section to be substantial 
- larger than typically provided by the effects of large-scale 
structure. So, one must be lucky to observe higher order 
imaging produced by such configurations. Six and eight im- 
age configurations can also occur when multiple galaxies lie 
at the same redshift (Keeton et al. 2000) or at distinct red- 
shifts (Chae, Mao & Augusto 2001). This may be quite com- 
mon, as the lens in the 6 image case found by Rusin et al 
(2001) seems to be a cluster. However, it is perhaps less sur- 
prising that groups of galaxies can produce unusual image 
configurations. 

A crude estimate suggests that about 3% of the lensing 
galaxy population has the potential to produce sextuplets 
or octuplets. Unfortunately, the cross-sectional area in the 
source plane producing sextuple imaging is approximately 
a factor of 10 times smaller than that producing quadru- 
ple. However, the probability of observing a sextuplet is in- 
creased by a factor of 5 by the magnification bias (see 
e.g. Turner 1980; Turner, Ostriker & Gott 1984). There- 
fore, a realistic estimate is that sextuple and octuple images 
comprise ~ 1% of all multiply imaged systems. Note that 
we have ignored galaxies that show isophote twisting (e.g., 
Bertola 1981) and galaxies which exhibit boxiness and disk- 
iness in the same isophotal contours (e.g., Nieto & Bender 
1989). We have also not included the contribution from mul- 
tiple lensing by groups of galaxies. These might be expected 
to augment the numbers of sextuplets and octuplets, and so 
we reckon our estimate errs on the side of pessimism. 

It is reasonable to expect that both the Next Genera- 
tion Space Telescope (NGST; see e.g. Barkana & Loeb 2000) 
and the Global Astrometric Interferometer for Astrophysics 
(GAIA; see e.g., de Boer et al. 2000) will discover such con- 
figurations. For example, GAIA is the ESA satellite now 
selected as a Cornerstone 6 mission as part of the Science 
Program ^. It is an all-sky survey satellite that provides 
multi-colour, multi-epoch photometry, astrometry and ra- 
dial velocities on all objects brighter then V ~ 21 (e.g., de 
Boer et al. 2000). It will be sensitive to multiply-imaged 
quasars with separations as small as ~ 0.2". It is crudely 
estimated that GAIA will detect ~ 4000 multiply-imaged 
quasars. In a dataset of such size, we expect ~ 40 sextuple 
or octuple image systems. Both NGST and GAIA will also 
find "mini" Einstein rings. We expect that, in many cases. 



t http:/ /astro. estec.esa.nl/SA-general/Projects/GAIA/ 
gaia.html 
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some of the six or eight images will be merged to form arcs. 
On account of the magnification bias, we expect the lenses to 
be typically high redshift galaxies and so the sources them- 
selves will have redshift > 2. Therefore, GAIA should also 
find "mini" Einstein rings, where the diameter of the ring is 



sub- arcbecuud 



It 13 intcrGsting to note that any additional cusps pro — 
duced by distortion or deviation of elliptical symmetry in- 
creases the magnification of the quadruple image configura- 
tions as well. It is still a puzzle why we observe so many 
quadruplets. Simulations typically predict many more dou- 
ble than quadruple image systems. Boxy or disky ellipti- 
cal galaxies favour the occurrence of more quadruplets This 
is because the magnification of the 4 images is increased 
by the appearance of additional cusps inside the 4 image 
regime. This fact has been neglected as only highly sym- 
metric galaxy models have been typically used in numerical 
simulations. 

The boxy and disky galaxy models with flat rotation 
curves introduced in this paper are a natural extension of the 
popular and widely-used elliptic density models. They have 
many attractive properties. The critical curves and the rela- 
tive time delays are analytic. The deflection angles and the 
caustics are simple quadratures and so easily evaluated. The 
lens equation can be reduced to a one-dimensional equation, 
and so the image positions can be found numerically with- 
out recourse to grid searching algorithms. Boxy and disky 
galaxy models provide new opportunities to remodel ob- 
served quadruple lenses, where existing models have proved 
unable to calculate the flux ratio of the images correctly. A 
good candidate might be MG0414-I-0534 (see e.g. Hewitt et 
al. 1992). Here, we have 4 radio images with well-determined 
flux ratios; however, no model has so far been able to repro- 
duce the flux ratios satisfactorily (Falco, Lehar & Shapiro 
1997). The most recent attempt (Trotter, Winn & Hewitt 
2000) concludes that the lensing galaxy is somewhat quad- 
rangular and boxy. 
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